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ABSTRACT 

The young pulsar PSR B0540-69 was one of the first targets observed with the Chandra 
X-Ray Observatory. The high angular resolution of Chandra allows us to resolve the compact 
nebula surrounding the pulsar. We have determined a position for PSR B0540-69 of R.A. = 
05''40"nl^221,decl. = -69°19'54".98 (J2000) with a la radial uncertainty of 0.7". Combining 
our measurements of the pulsar period with previous measurements covering a span of 12 years, 
we derive a braking index of 2.08. The spectrum of the pulsed emission is consistent with a 
power-law with a photon index of 1.83 ± 0.13. The compact nebula has a softer spectrum with 
a photon index of 1.85-2.26. 

Subject headings: pulsars: general — pulsars; individual (PSR 0540-69) — stars: neutron — supernova 
remnants — X-rays: stars 



1. Introduction 

PSR B0540-69 is a young, energetic pulsar lo- 
cated in the Large Magellanic Cloud (LMC) and 
discovered in soft x-rays using the Einstein ob- 
servatory (Seward, Harnden, & Helfand 1984). 
The pulsar is similar to the Crab pulsar in pe- 
riod, 50 ms versus 33 ms for the Crab, charac- 
teristic spin-down age, 1600 yr versus 1240 yr, 
and spin-down power, 1.5 x lO'^^ergs"^ versus 
4.7 X lO^^ergs"^. Comparison of detailed obser- 
vations of PSR B0540-69 versus those of the Crab 
should be a useful test of models of young pulsars, 
their evolution, and their associated nebulae. 

Due to the large distance to the LMC, PSR 
B0540-69 is a faint source. While it has been been 
observed in the x-ray, optical, and radio, relatively 
long observations on large telescopes are required 
to obtain good signal to noise and, thus, observa- 
tional coverage is spotty. Previous observations of 
PSR B0540-69 have led to discordant ephemerides 
and significant disagreements about the braking 
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index. 

Here, we present an analysis of early Chandra 
X-Ray Observatory observations of PSR B0540- 
69 made using both the High Resolution Camera 
(HRC; Murray et al. 1997) and the ACIS imag- 
ing camera (ACIS-I; Bautz et al. 1998). An early 
analysis of the HRC data has been previously de- 
scribed in Gotthelf & Wang (2000) which reported 
on the morphology of the spatially resolved neb- 
ula near the pulsar and outer shell of the remnant. 
Here, we analyze the HRC data with an improved 
aspect solution needed to obtain the correct po- 
sition for PSR B0540-69 and to remove artifacts 
present in the initial processing of the data. In 
addition, we present results on the x-ray spectrum 
of the pulsar, the compact nebula, and the outer 
remnant from an analysis of the ACIS data. We 
describe the observation and our analysis in §2. 
We present results on the source location, the pul- 
sar period history, the extent of the nebula, and 
the spectrum of the emission in §3. We conclude 
in §4. 

2. Observations and Analysis 

PSR B0540-69 was observed with the Chan- 
dra X-Ray Observatory (CXO; Weisskopf 1988) 
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during the verification and checkout phase on 31 
August 1999 using the High Resolution Camera 
(HRC; Murray et al. 1997) for a total of 18 ks 
of good observing time and on 26 August 1999 
using the ACIS imaging camera (AClS-1) for a to- 
tal of 20 ks of good observing time. All obser- 
vations employed the High- Resolution Mirror As- 
sembly (HRMA; van Speybroeck et al. 1997). 

The absolute time calibration of Chandra has 
not been finalized. For this reason, we use only fre- 
quency measurements in the analysis below. The 
accuracy of the relative time-tagging within an 
observation is determined by the spacecraft clock 
which is stable to better than one part in 10^ over 
one day (Chandra Observatory Proposer's Guide). 
The conversion to a barycentric time frame intro- 
duces an additional frequency uncertainty of less 
than 2 x 10~^ in these observations. Thus, the ac- 
curacy of the barycentric frcqucnc;y determination 
is more than sufficient for the frequency analysis 
presented below. 

2.1. HRC Analysis 

The HRC is a microchannel plate imager hav- 
ing very good spatial and time resolution, but es- 
sentially no energy resolution. Each photon de- 
tected by the HRC is time tagged with a precision 
of 16 /is and position tagged with a precision of 
0.132". This makes possible time-resolvcd imag- 
ing studies. After launch, an error was found in 
the HRC wiring which causes the event time tag 
to actually contain the time of the previous event 
trigger (Seward 2000). Without correction, this 
leads to an timing error equal to the time between 
successive event triggers. The event trigger rate 
is typically 250-300 s~^, leading to a typical er- 
ror of 3-4 ms. A special operating mode for the 
HRC has been developed in which all event trig- 
gers are telemetered to the ground and the timing 
error can be eliminated by shifting of event time 
tags. However, in the HRC observation reported 
here, not all event triggers result in valid events 
telemetered to the ground and correction of the 
event times is not possible in all cases. Thus, the 
event times contain typical errors of 3-4 ms. This 
error is small relative to the 50 ms period of PSR 
B0540 69 and does not significantly affect the pe- 
riod determination or the pulsed-phased resolved 
imaging presented below. 

This observation of PSR B0540-69 was only the 



third observation made with the HRC of a celes- 
tial source and was performed before the correct 
focus position of the HRC was determined. Dur- 
ing this observation, the HRC was displaced by 
260 /um from the true focus. This defocuses the 
image of an ideal point source, ignoring mirror and 
aspect imperfections, to an annulus with inner and 
outer radii of 0.17" and 0.33". 

We applied aspect to X-ray events from the 
HRC using a version of the Chandra X-Ray Obser- 
vatory Center (CXC) aspect pipeline (R4CU5, 11- 
Jan-2000) which utilizes elliptical Gaussian cen- 
troiding. Use of the improved aspect solution is 
critical to obtain correct positions. The relative 
aspect reconstruction within the observation was 
checked by using the aspect solution for this obser- 
vation to de-dither the optical guide star centroids. 
This is similar to the process used to apply as- 
pect correction to x-ray events. We found no sys- 
tematic residuals in the de-dithered star images, 
leading us to conclude that relative aspect during 
the observation is stable to within 0.3". The ab- 
solute accuracy of the aspect reconstruction was 
tested using a series of observations of 8 HRC- 
I observations of known x-ray point sources with 
identified optical counterparts with very accurate 
astrometry (Aldcroft et al. 2000). In each case, 
the derived X-ray source position was compared 
with the known counterpart position to derive the 
celestial location error. The absolute source coor- 
dinate uncertainty is 0.7" {la radial error). 

The data were filtered using event screen- 
ing techniques (Murray et al. 2000) to eliminate 
"ghost" events produced by the HRC electronics 
(Dobrzycki 2000). Screening to eliminate "ghost" 
events is particularly important due to the bright- 
ness of PSR B0540-69. The screening removes 
image artifacts which could have been incorrectly 
interpreted as features of the x-ray emission. 

We extracted an image which shows the pul- 
sar as a bright point source surrounded by nebu- 
lar emission with a roughly elliptical shape with 
an extent of a few arc-seconds. We used cen- 
troiding as implemented in the standard Chan- 
dra software routine celldetect to determine the 
position of the point source (CIAO VI. 1 Software 
Tools Manual). The coordinates found are R.A. = 
05M0™llM68,decL = -69°19'55".12 (J2000). 

These coordinates were used to perform an ini- 
tial barycentric correction employing the JPL DE- 
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405 cphcmcris. Wc searched for pulsations using 
the Z"^ statistic also known as the Rayleigh test 
(e.g. Buccheri, Ozel, & Sacco 1987). As the pulse 
profile of PSR B0540 69 is known to be roughly 
sinusoidal (Seward et al. 1984), we used only the 
first order Z'^ statistic. Pulsations were easily de- 
tected from PSR B0540-69. 

As the diffuse emission from the nebula sur- 
rounding PSR B0540-69 may affect the position 
determination, we created a difference image to 
isolate the pulsed source. We made a "pulsar- 
on" image using photons taken from an interval 
of 0.5 in phase with an offset chosen to maximize 
the number of photons selected, and a "pulsar-off " 
image from the remaining data. The "pulsar-off" 
image was subtracted from the "pulsar-on" image 
to obtain the difference image. All x-ray sources, 
in particular the constant nebular emission, not 
pulsed at the folding frequency are removed by 
this procedure. The difference image shows only a 
single point source. 

We note that the "pulsar-off" image does con- 
tain some pulsar emission because the pulse pro- 
file is roughly sinusoidal (Seward et al. 1984) and 
is not zero during the full 0.5 in phase of the 
"pulsar-off" image. The HRC timing error causes 
some additional pulsar photons to be shifted into 
the "pulsar-off" image, but the effect is relatively 
small since the timing error of 3-4 ms is small com- 
pared to the phase bin size of 25 ms. The pres- 
ence of some pulsar emission in the "pulsar-off" 
image reduces the amplitude of the pulsar peak in 
the difference image. However, we found that the 
difference image still contains a highly significant 
peak at the pidsar position and that use of the 
maximum size phase bins was advantageous in re- 
ducing image noise away from the peak where the 
photon statistics arc much lower. We also note 
that the timing error produces a very small net 
shift of nebular photons into the "pulsar-on" im- 
age. If the timing error were simply a random 
number added to the time tag, then the net shift 
would be zero since the number of nebular pho- 
tons moved from the "pulsar-off" image into the 
"pulsar-on" image would equal the number moved 
in the opposite direction. However, since the er- 
ror is actually a shift of time tags, the fact that 
the event trigger rate is slightly higher during the 
"pulsar-on" interval (by a factor equal to the pul- 
sar count rate divided by average total event trig- 
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Fig. 1. — Contour plots of emission from PSR 
B0540-69. The solid contours are the "pulsar-off" 
phase interval and the dotted are the difference 
image. The cross is the pulsar position. Both im- 
ages were smoothed with a Gaussian with FWHM 
= 0.7". The contour levels are 0.1, 0.3, 0.5, 0.7, 
and 0.9 of the peak pixel. The dashed line and 
arrow labeled "cc" indicate the orientation of the 
ACIS image. 

ger rate or roughly 0.3%) produces a net shift. We 
estimate that the net shift of nebular photons into 
the "pulsar-on" image is less than 0.05%. 

Using the position obtained from this difference 
image, we performed a second barycentric correc- 
tion and selected events within a circle with a ra- 
dius of 9 HRC pixels (1.1") around the difference 
image position. We again searched for pulsations 
and found a strongly significant pulsation signal at 
19.7988001(21) Hz at epoch MDJ 51421.6240; the 
digits in parentheses indicate the \a uncertainty. 

We made a second difference image using this 
best period, see the dotted contours in Fig. 1. The 
position of the source in the two difference images 
agree within 0.3 pixel = 0.04". The position from 
the final image is R.A. = 05'^40'"11^221, deck = 
-69°19'54".98 (J2000). The position uncertainty 
is dominated by the accuracy of the aspect recon- 
struction which we take to be 0.7" (la radial er- 
ror). For the difference image, 50% of the flux 
is contained with a diameter of 1.3". This is sig- 
nificantly larger than the half-power diameter of 
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0.76" measured for AR Lac in calibration observa- 
tions made to determine the on-orbit point spread 
function (Jerius et al. 2000), but is consistent with 
the expected combined effects of the HRMA point- 
spread function given the harder spectrum of PSR 
B0540-69, the relative aspect accuracy measured 
for this observation, the position resolution of the 
HRC, and the defocusing. 

2.2. ACIS Analysis 

The ACIS-I camera consists of an array of 4 
front-illuminated charged coupled devices (CCDs). 
The physical pixel size is 0.24 /xm which at the 
aim-point of the telescope is comparable to the 
0.5" spatial resolving power of the HRMA. Each 
CCD contains 1024 x 1024 pixels organized into 4 
readout nodes each of which reads out 1024 rows 
and 256 columns of pixels. The pulsar was posi- 
tioned at the aimpoint of the telescope, which is 
about 960 rows away from the readout node on 
ACIS-I device 13. The pulsar and its surround- 
ing nebula fit very well onto node 3 of device 13, 
therefore we use only this node in our analysis. 
The ACIS flight software was set to an on-board 
processing mode that retains the 3x3 pixel event 
island. 

In order to have sufficient time resolution to re- 
solve the pulses from PSR B0540-69, the ACIS-I 
was operated in continuous clocking mode ("cc- 
mode" ). In this mode, charge is continually shifted 
between successive pixels in each row and the pixel 
at the base of each row is read out after each 
charge transfer. This leads to acquisition of a one- 
dimensional image each 3.1 ms. The 1-d image 
integrates flux along the readout direction. The 
time of arrival of each photon must be calculated 
from the time of readout and the (unknown) posi- 
tion of arrival along the readout direction divided 
by the speed of pixel transfer along the readout 
direction. 

The data on ground were manually processed 
through the standard CXC pipelines into CXC 
level 1 event lists. As software for the analysis of 
cc-mode data is not currently available from the 
CXC, we further processed the ACIS-I data us- 
ing custom IDL procedures. The pointing direc- 
tion was transformed to the spacecraft coordinate 
frame {sy, s^, in detector pixels) so that the event 
positions could be corrected for dither in the detec- 
tor X direction only via Xp = Xt + S * 5x{ti) + Xr, 




6 4 2 -2 -4 -6 
Arcsec from pulsar 




40 20 -20 -40 
Arcsec from Pulsar 



Fig. 2. — One-dimensional images from ACIS-I 

cc-mode data. The top panel shows two image 
profiles. The solid curve shows events with pulse 
phases between 0.1 and 0.3 where the pulse is max- 
imiim and the dashed curve shows phases 0.6 to 
0.8 where the pulse is at minimum. The bottom 
panel shows the overall profile. The dotted lines 
indicate the various regions used in the spectral 
analysis. 

where Xp is the position of the event projected 
along a line on the sky (at a position angle of 
66.97°), Xt is the event coordinate in "tiled" co- 
ordinates (tdetx), S is the sign of the correction 
for dither (-1 for chips II and 13 and -1-1 for chips 
10 and 12), and Xr is a random variable between 
-0.5 and 0.5 that reduces effects of aliasing. The 
spacecraft dither and optical bench distortions are 
accounted in the term Sx{t), which is the derived 
from interpolating Sx{t) — Sz{t) — dzsiM{t), where 
dzsiM{t) gives the motion of the science instru- 
ment module (SIM) measured by the on-board 
fiducial lights. 

The histogram of the Xp values, i.e. the one- 
dimensional image perpendicular to the readout 
direction, is shown in Fig. 2. For this observa- 
tion, the 1-d image is along a positive angle of 
67°, or approximately ENE, as indicated on the 
HRC image in Fig. 1. We define the angle along 
the 1-d image as positive to the ENE. A reference 
point is computed by fitting a narrow Gaussian 
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to the brightest pixels, assuming that the syn- 
chrotron nebula provides a simply sloping back- 
ground. The profile of the pulsar is well fitted by 
a Gaussian with a dispersion of 0.30", which is 
consistent with the telescope point spread func- 
tion (PSF). The coordinates in Fig. 2 are relative 
to the fitted peak position. The sharp peak of the 
pulsar is flanked by the integrated image of the 
compact nebula. The compact nebula brightness 
distribution peaks to the SW from the pulsar po- 
sition and declines rapidly further SW after the 
peak. The brightness falls off more slowly toward 
the NE. The compact nebula gives way to a more 
diffuse remnant, the outer supernova remnant, at 
about 5" from the pulsar. The diffuse emission 
extends to ~ 35" in one direction and is truncated 
by a chip gap at 50" from the pulsar in the other 
direction. 

The event times were constructed using the 
event row numbers, t/,, an interpolation of the 
dither and SIM correction, Sy{t) = Sy{t) — 
dysiM{t); the ACIS exposure time to UT lookup 
table, fc/T(e); the geocenter and solar system 
barycenter corrections, dt^ and dt^; and a ran- 
dom variable, yr, between and 1 that reduces 
effects of aliasing: tTT,i = tuT{e-i) + tfiHi + 
Vr) + {tTT — tur), where the last term is sim- 
ply a constant for this observation (64.184 s), 
tf is the average detector frame shift time in 
ephemeris time units (about 0.00284996 s). The 
exposure time to UT lookup table should be ac- 
curate to less than 1 msec. Finally, the times 
corrected to the solar system barycenter (TCB) 
are U = txTA ~~ Sy{tTT,t) +dt(i,{tTT,i) +dtQitTT,i)- 

We have found that an additional 5 s must be 
added to tuT,i in order to match the TT times 
computed for the aspect data. The 5 s offset is 
a data processing anomaly that has not been ex- 
plained fully so we merely report that it exists. 
Specifically, the offset is needed to bring the time- 
tagged aspect data into concurrence with the time- 
tagged photon events. Since we believe that the 
aspect data have "correct" times (based on data 
processing of other instrument configurations) , we 
have adjusted the event times. We can state that 
this 5 s offset is good to about 1 s and that it has 
not changed across the observation at this level. 
More importantly, different cc-mode observations, 
even ~ 50 days later, give the same value of the 
shift to within 1 s, thus the secular variation in 



the offset is no larger than 3 x 10~ . If there were 
a secular drift of this timing offset at this level, 
then it would produce a 5 ms timing error across 
the ACIS observations and lead to an error on the 
period roughly twice the statistical error quoted 
below. We use this larger error in the ephemeris 
fit described below. The barycenter timing correc- 
tion changes by only 60 ms across the ACIS ob- 
servation and is approximately quadratic in form, 
thus a 5 s timing uncertainty produces an uncer- 
tainty of less than 1.4 x 10~^ on the barycenter 
frequency correction. 

The events within 0.7" of the peak in the 
Xp distribution were Fourier transformed to ob- 
tain the candidate pulse frequency and the best 
estimate was determined from pulse folding: 
19.7988900(30) Hz. Two imaging profiles were 
obtained by selecting from phases 0.6 to 0.8 where 
the pulse is at a fiat minimum, and phases 0.1 to 
0.3 where the pulse is nearly level at maximum. 
These two image profiles, along with the profile of 
all the emission, are shown in Fig. 2. 

2.2.1. Spectral analysis 

The observation took place in a period where 
the ACIS instrument experienced large changes in 
the spectral response due to low energy cosmic 
proton impacts. This resulted in a significantly 
increased charge transfer inefficiency (CTI). This 
has two effects on the spectral response. First, 
the energy scale changes due to charge loss that 
is not recovered during the readout, and, second, 
the introduction of additional electron noise due 
to enhanced charge trapping worsens the spectral 
resolving power. 

The focal plane temperature of the device at 
the time of the observation was -100° C. Unfortu- 
nately, during this period of change there are no 
response matrices available at a focal plane tem- 
perature of — 100°C. However, there are calibra- 
tion products available for a focal plane temper- 
ature of — 110°C. Calibration data indicate that 
throughout most of the energy scale the change in 
CTI with focal plane temperature is fairly linear. 
The change in spectral resolution is then assumed 
to be linear as well. Calibration data also indi- 
cate how CTI changed during the period of degra- 
dation. By estimating the change in energy scale 
due to CTI as well as taking in account the energy 
scale shift implied by the temperature change from 
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Fig. 3.— Spectrum from PSR B0540-69. 

-100° C to -110° C, we find that the detector area 
between row 100 and 320 at —110° C matches the 
energy scale of the detector node at the time of the 
observation between rows 800 and 1024 at —100° C 
quite well. In order to fine tune the energy scale, 
we made fits using the available response matrices 
at -110° C valid for rows 100 to 320 and selected 
the ones that matched the expected mirror Ir-odgc 
structures and position with the one expected for 
the model effective area. Currently available re- 
sponse matrices show a quite strong non-linearity 
in the energy scale below 1 keV that can affect 
the analysis of the soft part of the spectrum. For 
the analysis here, wo do not fit the spectra be- 
low 1.2 keV. Fig. 3 shows a fit of the pulsar spec- 
trum (see below). By examining known edges in 
the CCD response, particularly the Si K edge at 
1.738 keV, we estimate that the energy scale is 
accurate to within 2%. 

Background analysis is important in cc-mode 
because each pixel read out has the integrated 
background from 1024 CCD pixels. The contri- 
bution does not exceed 2% for the pulsar or its 
compact nebula, but is almost 30% for the outer 
remnant. We use the source free region between 
—35" and —50" to compute the backgromid spec- 
trum after cleaning the event list for all apparent 
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Fig. 4. — Deviation of frequency from a linear 
ephemeris. The crosses indicate data from Deeter 
et al. (1999) and the solid line is the correspond- 
ing quadratic ephemeris. The four diamonds in- 
dicate more recent data from, in chronological or- 
der, Boyd et al. (1995), Mineo et al. (1999), and 
the ACIS and HRC values from this work. 

hot columns. The measured average bac;kground 
rate was 1.20 ± 0.01 x 10"-'' ctss^^ cohmm"^ 

For the effective area we used the one avail- 
able for the standard ACIS grade set, which is 
appropriate since the observation was performed 
in a mode that allowed reconstruction of a 3 x 3 
pixel island for each event. In order to compute 
the corresponding ancillary data file (arf) we used 
custom software that averaged over the aspect so- 
lution. This is legitimate for our analysis because 
no source portion in the analysis crosses a CCD or 
node boundary. The spectral fits were performed 
with XSPEC v.10.0. 

3. Results and Discussion 
3.1. Timing 

In Fig. 4, we compare our two period measure- 
ments and two other recent observations of PSR 
B0540-69 (Boyd et al. 1995; Mineo et al. 1999) 
with the ephemeris derived by Deeter et al. (1999) 
using data from MJD 46992-48335. Extrapola- 
tion of the Deeter et al. (1999) ephemeris shows 
that the new data, which extend almost 8 years 
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Table 1 

Spin Parameters for PSR 0540-69. 



Parameter 


Value 


to (A'lJD) 


47700.0 


z^o (Hz) 


19.8593584982(40) 


z>o (10-i"Hzs~i) 


-1.8894081(7) 


i>o (10-21 Hz 


3.7425(43) 


n (braking index) 


2.0820(24) 



Note. — Numbers in parentheses are la 
errors in the last quoted digits. 



beyond the data used to derive the ephemeris, dif- 
fer by at most 2a from the predicted frequencies. 
We interpret this as a strong indication that the 
timing solution found by Deeter et al. (1999) is 
accurate. Both Chandra frequencies are consis- 
tent, within the la errors, with the extrapolation 
of the Deeter et al. (1999) ephemeris. The HRC 
frequency measurement is within 1.9 x 10~^Hz of 
the predicted value and suggests that no large per- 
sistent offsets in pulse frequency or its first deriva- 
tive occurred during the 8 year interval. However, 
we are unable to rule out small persistent offsets 
Az/ < 2 X 10"^ Hz or AzX 8 X IQ-i^Hzs-i, sim- 
ilar to those reported by Deeter et al. (1999). 

We combined our frequency measurements with 
those of Boyd et al. (1995), Minco et al. (1999), 
and Deeter et al. (1999; specifically the frequency 
data Tiscd to derive the ephemeris and marked 
with note 5 in Table 4) to derive the ephemeris 
in Table 1. The ephemeris is consistent, within 
errors, with that of Deeter et al. (1999). We find 
a braking index of n = 2.0820 ± 0.0024 {la uncer- 
tainty) at an epoch 47700 M.JD. This braking in- 
dex is much lower than the index of 3 expected for 
energy loss via electromagnetic dipole radiation 
and is significantly lower than the indices mea- 
sured for the Crab pulsar or PSR 1509-58. The 
low value of the braking index may be due to a 
pulsar wind (Manchester & Peterson 1989; Na- 
gase et al. 1990). The synchrotron nebula may 
provide information on the properties of the pul- 
sar wind and it would be of interest if the prop- 
erties of nebulae surrounding young pulsars could 



be related to their braking indices. Other possible 
reasons for the low value of the braking index are 
distortion of the magnetic field lines (Manchester 
& Taylor 1977), or a time- varying magnetic field 
strength (Blandford & Romani 1988). Continued 
pulse timing of PSR B0540 69 should make it pos- 
sible to distinguish amongst these alternatives. 

3.2. Position 

Positions for PSR B0540-69, or its optical 
counterpart, have been reported by Seward et 
al. (1984) from X-ray imaging, by Deeter et al. 
(1999) from X-ray and optical timing, by Caraveo 
et al. (1992) and Schmidtke et al. (1999) from op- 
tical imaging, and by Shearer et al. (1994) from 
time-resolved optical imaging. Our position is in 
good agreement with those of Schmidtke et al. 
(1999) and Shearer et al. (1994) and marginally 
consistent with that of Caraveo et al. (1992). The 
results of Shearer et al. (1994) indicate that the 
pulsar may be slightly displaced from the maxi- 
mum in the steady optical emission identified with 
the pulsar by Caraveo et al. (1992) in a direction 
which would improve the agreement with our po- 
sition measurement. Our HRC position for PSR 
B0540-69 lies 2.6" from the original X-ray posi- 
tion reported by Seward et al. (1984), although 
within the 90% confidence error box. Our posi- 
tion is 1.8" from the position derived from X-ray 
timing measurements with the Ginga satellite by 
Deeter et al. (1999). Interestingly, the Deeter et 
al. (1999) position agrees very well with the other 
positions determinations in declination, but is sig- 
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nificantly displaced from all the other positions in 
R.A. The location of PSR B0540-69 less than 4° 
from the ecliptic pole may lessen the accuracy of 
the determination of its R.A. from pulse timing. 

Because a shift of 2.6" in source position will in- 
duce an annual sinusoidal variation in the barycen- 
tric pulse arrival times with an amplitude of sev- 
eral milliseconds, an accurate position is critical 
when performing a pulse time of arrival (TOA) 
analysis. Most timing studies of PSR B0540-69 
have used the original X-ray position and should 
be reconsidered in light of the more accurate po- 
sition reported here. In particular, the previous 
position uncertainty may have contributed to the 
disparate range of braking indices found for PSR 
B0540-69. Also, the decreased uncertainty in the 
position may improve the evidence presented for 
the detection of timing noise in PSR B0540-69 
(Eikenberry, Fazio, & Ransom 1998). 

3.3. Nebula 

We found a 0.35" displacement between the 
controid of the total x-ray omission near PSR 
B0540-69 and the pulsar itself. This relative dis- 
placement is significant and represents a physical 
offset between the pulsar and the peak x-ray emis- 
sion in the nebula similar to the offset seen in the 
Crab nebula versus pulsar. 

The difference image, which isolates the pre- 
sumably point-like pulsar emission, provides an 
indication of the point-spread function for a point 
source including the finite resolution of the HRMA 
and HRC, the effects of the aspect errors, and the 
focusing error described above. Comparison of the 
difference image and the "pulsar-off" image, which 
highlights the nebular emission by reducing - but 
not fully eliminating - the pulsar emission, shows 
that the nebula is spatially extended and has a 
centroid which is slightly WSW from the pulsar, 
see Fig. 1 in which the dotted contours represent 
the difference image and the solid contours repre- 
sent the "pulsar-off" image. 

To estimate the true extent of the nebula, we 
performed an image deconvolution using the differ- 
ence image described above to represent the tele- 
scope point-spread function including the effects of 
aspect errors and defocusing. We deconvolved an 
HRC image binned to 0.26" resolution using the 
standard Richardson-Lucy technique (Richardson 
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Fig. 5. — Deconvolved image of the emission from 
PSR B0540-69. North is up and the axes indicate 
displacement from the central pixel in arcseconds. 
The pixel size is 0.26". The scale indicates counts 
per pixel. The central pixel is saturated with 3230 
counts. 

1972; Lucy 1974). Application of the deconvolu- 
tion and also the event screening described above 
is essential to determine the true angular extent 
of the nebula. The resultant deconvolved image 
in shown in Fig. 5. The nebula has an angular 
diameter of 2" — 3" and is more extended along a 
NE-SW axis. 

Gotthelf & Wang found the same orientation for 
the compact nebula, but a larger size. Their re- 
sults arc consistent with our un-dcconvolved con- 
tour plot in Fig. 1. The deconvolution gives a bet- 
ter estimate of the true angular extent of the neb- 
ula. Gotthelf & Wang (2000) also report marginal 
evidence for a "jet-like feature" extending to the 
NW from the pulsar. A similar weak feature 
is present in our data. However, the feature is 
present in the "pulsar-on" and difference images, 
but is significantly weaker in the "pulsar-off" im- 
age. Caution is required before interpreting this as 
a physical feature and another observation, with 
a different roll-angle, would be required for confir- 
mation. 
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The total HRC count rate from the nebula and 
pulsar within a 4" radius is 0.84 counts s~^. As- 
suming a power-law spectral model with a photon 
spectral index of 2.0 and an interstellar absorption 
column density Nji = 4.6 x 10^^ cm~^ (see the 
spectral fitting below) gives an unabsorbed flux of 
5.8 X 10"^^ ergcm-2 s"^ in the 0.2-10 keV band. 

3.4. Pulsed emission 

The HRC count rate within 1.1" of the pulsar is 
0.38 counts with 41% of the signal pulsed. For 
a narrower radius of 0.6", the pulsed fraction in- 
creases to 55%; however, some of the pulsed emis- 
sion is lost due to the large PSF during the HRC 
observation. We note that the measured pulse 
fraction is sensitive to the shape of the pulse near 
the maximum and minimum, which may be af- 
fected by the 3-4 ms timing error in the HRC. 
Thus, the true pulsed fraction may be somewhat 
higher than this value. 

This pulsed fraction is significantly higher than 
previously reported, e.g. ~ 15% from ROSAT 
(Finley et al. 1993), because Chandra can resolve 
the surrounding compact nebula. Using a power- 
law spectral model with a photon spectral index of 
1.83 and a column density Nh ~ 4.6 x lO^'^ cm~^ 
(see the spectral fitting below) gives an unab- 
sorbed pulsed flux of 1.1 X 10~^^ ergcm"^ s~^ in 
the 0.2-10 keV band. The off-phase deconvolu- 
tion (not shown here) shows only a weak point 
source at the pulsar position, indicating that, like 
the Crab Pulsar, any steady emission from the pul- 
sar itself is small compared with the pulsed com- 
ponent. 

The ACIS pulse profile for photons within 0.7" 
of the pulsar position is shown in Fig. 6. The 
ACIS cc-mode data suffer from a significant neb- 
ular contamination of the pulsar signal due to the 
integration of counts along the readout direction. 
Using the analysis of the position profile described 
above, we estimate the nebular contribution to 
be 570 counts per phase bin, leading to a pulse 
fraction of 62 ± 5%. This estimate is consistent 
with that found from the HRC data. Fig. 6 also 
shows the pulse profile for photons with energies 
above 4 keV and below 1.5 keV. The modulation 
increases with energy indicating that the pulsed 
emission is harder than the unpulsed or nebular 
emission. The pulsed fraction, even below 1.5 keV, 
is significantly higher than reported from Rosat 




0.4 0.6 
Pulse phase 



Fig. 6. — Pulse phase histogram of the events from 
PSR B0540~69 in the ACIS-I data. The solid line 
indicates photons within 0.7" of the pulsar posi- 
tion. The dashed line is for photons with energies 
above 4 keV and the dotted line is for photons 
with energies below 1.5 keV; both normalized to 
have the same total number of counts as in the full 
band. Representative error bars are shown. 



(Finley et al. 1993) or BeppoSAX (Mineo et al. 
1999) due to nebular contamination in the previ- 
ous results. 

3.5. Spectra 

We accumulated ACIS spectra from the pul- 
sar, two regions surrounding the pulsar where the 
compact nebula is prominent, and the outer re- 
gion where the outer remnant dominates. These 
spectra were analyzed using the response matrices 
described above. In this section, we quote x-ray 
fiuxes for the energy band 0.6-10 keV and errors 
for 90% confidence unless otherwise noted. 

3.5.1. Pulsar spectrum 

For the pulsar emission we selected events 
within 0.7" of the peak position in Fig. 2. The 
background subtracted source rate was 1.49 cts/s. 
Pile-up here is negligible because of the small 
clocking time in cc-modc. We used this overall 
spectrum to determine the proper detector re- 
sponse matrix for the analysis. Since we ignore 
energies below 1.2 keV, we are less able to con- 
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Table 2 
Parameters of Spectral Fits. 





Index 


Flux 




Pulsar all 


1.92±0.11 


1.78±0.05 


0.99 


Pulsar "off" 


2.09±0.14 


0.21±0.08 


0.60 


Pulsar "on" 


1.88±0.11 


0.74±0.04 


0.75 


Pulsed 


1.83±0.13 


0.53±0.06 


0.75 


Nebula (left) 


1.96±0.11 


0.98±0.01 


1.14 


Nebula (right) 


2.12±0.14 


0.49±0.01 


0.63 



Note. — This Table contains the photon in- 
dex from the power- law spectral fit, the ab- 
sorbed flux in the 0.6-10 keV band in units of 
10^^^ ergcm^^ s~^, and the xt of the spectral fit. 
The errors are 90% confidence. 



strain the column density, Nh, in our fits. In 
order to fit all spectra equally well, it was best to 
fix Nh = 4.6 X lO^^cm"^. The overall spectrum 
(see Table 2) is well represented by a power law of 
index 1.92 ±0.11. Both Nh and the photon index 
are in good agreement with the values found by 
Finley et al. (1993) with the Rosat PSPC and by 
Mineo et al. (1999) with BeppoSAX. 

In order to measure the spectrum of the pulsed 
emission alone, we have to subtract the back- 
ground imposed by the emission of the nebula. 
This can be done by making use of the pulse phase 
resolved components. We accumulated spectra 
at pulse phases 0.6-0.8 when the pulsed emis- 
sion is at a minimum, the "pulsar-off" spectrum, 
and at pulse phases 0.1 0.3 when the pulsed emis- 
sion is maximum, the "pulsar-on" spectrum. We 
subtracted the "pulsar-off" spectrum from the 
"pulsar-on" spectrum to produce the pulsed emis- 
sion spectrum. Table 2 shows the result of the 
spectral fits. The "pulsar-off" spectrum has a 
somewhat softer spectrum than the "pulsar-on" 
spectrum. The pulsed emission spectrum shows 
an even harder spectral index of 1.83 ± 0.13. The 
total absorbed pulsed flux can be calculated by 
integrating the emission above the lowest pulse 
phase. If we consider the flux level at phase 0.7 
as the unpulsed flux level, we find a pulsed flux of 



7.4 ± 0.2 X 10"-^^ ergcm"^ s~^. This corresponds 
to an unabsorbed flux of 1.2 x 10~^^ ergcm"^ 
in the 0.2-10 keV band, and is in good agreement 
with the HRC result above. Both the pulsed flux 
and the photon index of the pulsed emission are in 
good agreement with previous results from Rosat 
(Finley et al. 1993) and BeppoSAX (Mineo et al. 
(1999). 

3.5.2. Nebula spectrum 

Fig. 2 shows the pulsar flanked by the emis- 
sion from the nebula. We accumulated spectra on 
both sides of the pulsar, from 0.7" to 4.7" (left) 
and -0.7" to -4.3" (right) and fitted both spec- 
tra with a power law. The flux from the nebula is 
similar to that from the pulsar and thus the spec- 
tral parameters are similarly constrained. Here 
the fitted column density is somewhat lower, which 
again is a systematic effect of the low energy re- 
sponse below 1 keV, which now covers a somewhat 
larger region on the device. The spectral index of 
the nebula, when fixing Nh to the value used for 
the pulsar, is slightly steeper than that of the pul- 
sar, see Table 2. The total absorbed flux in the 
two flanks amounts to 1.47 x 10~^^ ergcm"^ s~^. 
The pulsar flux at low phase may also be asso- 
ciated with the nebula, which for the whole cy- 
cle amounts to 1.04 x 10~^^ ergcm"^ s~^. Thus, 
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Table 3 

Properties of the Crab and 0540-69. 





Crab 


0540-69 


Pulsar E (erg s ^) 


4.7 X 1038 


1.5 X 10^8 


Total Lx 


3.4 X 103^ 


2.1 X lO^'' 


Compact nebula Lx 


3.3 X 10^7 


1.7 X 10^7 


Pulsed Lx 


1.3 X 10^6 


0.4 X 10^7 


Nebula size (pc) 


1.5 X 1.5 


0.6 X 0.9 


Pulsar oflFset (pc) 


0.25 


0.09 


Braking index 


2.51 


2.08 



Note. — X-ray luminosities, Lx, are quoted in 
erg s~^ for the 0.2-10 keV band. 



the total absorbed flux for the compact nebula 
amounts to 2.51 ± 0.05 x 10"-^-^ erg cm~^ s~^. The 
corresponding total nebular unabsorbed flux is 
4.8 ± 0.05 x 10-" erg cm-^ s"! in the 0.2-10 keV 
band. This is good agreement with the unab- 
sorbed flux inferred from the HRC measurement. 

3.5.3. Outer remnant spectrum 

Figure 2 shows that the pulsar and it's com- 
pact nebula sits on top of some very extended ex- 
cess emission, the outer supernova remnant (Se- 
ward & Harnden 1994). We accumulated a spec- 
trum of this excess emission, in the region be- 
tween —4.8" and —30", and it turned out to be 
extremely soft, i.e. the bulk of the emission is 
found to be below 2 keV at a flux of approximately 
9.3 x 10^^3 erg cm"^ s~^. The current status of the 
ACIS response of front illuminated devices does 
not allow us to do spectral fits in this area yet, 
but there are line features evident in the range be- 
tween 0.7 and 1.8 keV and an overall resemblance 
to the spectrum observed in E0102-72 (Hayashi et 
al. 1994; Gaetz et al. 2000). 

4. Conclusion 

Study of PSR B0540-69 with its compact neb- 
ula, and comparison of its properties with those of 

the Crab and other young pulsars, should provide 
insights into the physical mechanisms for energy 



release in spin-powered pulsars including the dy- 
namics of the pulsar outflow. In Table 3, we com- 
pare the properties of PSR B0540-69 with those of 
the Crab, assuming a distance to PSR B0540-69 
of 55 kpc. 

The fraction of the total spin-down power in x- 
ray emission from the compact nebula and pulsar 
is 14% for PSR B0540 69 and half that for the 
Crab. The pulsed fraction of the total (pulsar plus 
compact nebula) flux is 20% for PSR B0540-69 
compared to only 4% for the Crab. In addition, 
the pulse shapes are quite different with the Crab 
showing a narrow pulse while PSR B0540-69 has 
a broad, almost sinusoidal, pulse. 

PSR B0540-69 has lower spin-down power and 
the surrounding nebula is less luminous. As a gen- 
eral rule (Seward & Wang 1985), the nebular lu- 
minosity varies with the pulsar spin-down power, 

• 1.4 

E, approximately as E ; the Chandra data are 
in rough agreement. The nebular volume for PSR 
B0540-69 is also considerably less than that of the 
Crab, although it is difficult to precisely deflne a 
"size" (the average FWHM is 0.5 pc and the av- 
erage FW at backgromid is 0.7 pc. The FW at 
background is given in Table 3). Since the nebula 
surrounding the pulsar is most likely due to syn- 
chrotron emission and energized by the pulsar, the 
smaller nebula likely reflects the lower spin-down 
power. However, the properties of the nebula may 
also depend on the surrounding environment and 
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the past history of the pulsar. 

The properties of the surrounding environment 
and also of the progenitor may be important in 
determining the nature of the outer remnant. In 
PSR B0540-69, we detect thermal x-ray emission 
from a surrounding remnant. This is in strong con- 
trast with the Crab pulsar, where no such remnant 
exists. The spectrum of the outer remnant of PSR 
B0540-69 appears to be line driven and thus more 
of the category we see in young oxygen-rich super- 
nova remnants of the Magellanic Clouds (Blair et 
al. 1999, Gaetz et al. 2000). 

PSR B0540-69 and its surrounding nebula are 
a remarkable parallel with the Crab Pulsar and 
its environment. It is almost certain that the 
same physical processes operate in both systems. 
However, further comparison shows several dis- 
tinct differences between the two systems, as re- 
vealed by the observations described above. Ad- 
ditional Chandra observations would allow us to 
study PSR B0540-69 in detail; making it possible 
to search for spectral structure in the synchrotron 
nebula, as has recently been detected from the 
Crab (Weisskopf et al. 2000). Confrontation of 
detailed models of the pulsar x-ray emission mech- 
anism, pulsar wind, synchrotron nebula, and the 
outer supernova remnant with these observations 
should help advance our understanding of spin- 
powered pulsars, their outflows, and the nature of 
their progenitors. 

We thank Rob Cameron for providing infor- 
mation on the orbit determination. We grate- 
fully acknowledge the efforts of the Chandra team 
and support from NASA Chandra contract NAS8- 
39073 and SAO grant SVl-61010. PK acknowl- 
edges partial support from NASA grant NAG5- 
7405. 
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